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Abstract 

The one-particle exchange in the double folded model is analyzed. To this aim the 
Extended Thomas-Fermi approach to the one-body density matrix is used. The nucleon- 
nucleon force with Yukawa, Gauss and Coulomb-type form factors are considered. The 
energy dependence of the exchange part of the double folded potential is investigated and a 
comparison of the present approach with former ones is carried out. 
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1 Introduction 

Up to the present moment the Double Folded Model (DFM) y], ||, (§], which starts from the 
effective nucleon-nucleon force and the single particle densities of the colliding nuclei, has become 
one of the most popular methods to calculate the real part of microscopic the Heavy-Ion (HI) 
optical potential. It is also known that folding methods are widely used in other problems of HI 
Physics H]]. While in the original DFM the antisymmetrization effects were taken into account 
by means of the effective zero-range pseudopotential , in the latest versions of the DFM the 
antisymmetrization has been explicitly considered ||, [3|, ||. It is assumed that the structure of 
of each isolated nucleus is described at Hartree-Fock (HF) level and that the exchange effects 
come from the non-diagonal part of the single-particle density matrix (DM). The intrinsic states 
of isolated nuclei are modified during the interaction. However in the DFM the so called "frozen 
density approximation" , which implies that the intrinsic states of the nuclei do not change during 
the interaction, is assumed. 

Thus the key point in the DFM is to express the non-diagonal part of the DM through its 
diagonal part or, in other words, in a Local Density Approximation (LDA) for DM. The simplest 
LDA is the Slater (SL) approach to the DM, which is equivalent to the Thomas-Fermi approach, 
valid in nuclear matter, and does not take into account the nuclear surface effects. A more 
elaborated LDA is given by the Density Matrix Expansion (DME) of Negele and Vautherin (NV) 
||. Another approach widely used in DFM [Q, |j is the Campi-Bouassy (CB) |?j approximation, 
in which by a proper choice of the effective momenta in the DME of NV, reduces to the SL form. 

The physical quantities that enter into the NV and CB approaches are just the Fermi mo- 
mentum fcp and the kinetic energy density r. In the HF method the single-particle density p as 
well as kF and r are determined selfconsistently. However, the attractive feature of DFM is to 
consider p as an input data. While the densities p can be obtained from the electron scattering 
data, kp and r should be defined theoretically. At present, the values of kp and r that corre- 
spond to a given single-particle density p are unknown. To determine them via p, approximate 
schemes are in order. In particular the Extended Thomas-Fermi (ETF) approach is widely used 
to this end: 




(1) 
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and 

5 2 3d p 3 

if a degeneracy 4 is assumed. 

Very recently the ETF approximation has been used for deriving the semiclassical one-body 
density matrix up to fi 2 order in the case of a non-local single-particle potential ||. This 
ETF DM is obtained starting from the Wigner-Kirkwood distribution function for a non-local 
one-body Hamiltonian ||. As is shown in ||, the NV and CB approaches to the DM with r 
calculated in the ETF approach, are just truncations of the full ETF DM. Consequently, it seems 
appealing to use the complete semiclassical DM to obtain the real part of the ion-ion potential 
and compare these results with former potentials obtained using these semiclassical NV and CB 
approximations. 

In the present paper we apply the ETF DM |8| for calculating the DFM potential. We obtain 
analytical expressions for the potential using the Gogny [jlOfl , M3Y [|J and Coulomb forces. To 
check the validity of our ETF approximation, we analyze the ie O- 16 reaction. The ground 
state of each isolated 16 nucleus is described using harmonic oscillator (HO) wavefunctions. 
In this case the quantal DFM potential can be calculated explicitly [111] and we compare this 
quantal potential with the semiclassical DFM potential obtained using the ETF, SL and CB 
approximations to the DM. 

The paper is organized as follows. In the first section we present the analytical derivation 
of the semiclassical DFM potential using the SL, CB and ETF approximations to the DM and 
Gauss, Yukawa and Coulomb form factors for the nucleon-nucleon force. In the next section 
we compare these semiclassical DFM nucleus-nucleus potentials obtained using several prescrip- 
tions for the density of eachisolated nucleus with the corresponding quantal potential. This 
comparison is done for the 16 0- 16 reactions. The summary and conclusions are given in the 
last section. 
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2 Formalism 



In order to obtain the HI potential in the DFM, one can start from its microscopic definition 

v(d) = [< »iM*i > - < >]> ( 3 ) 

iei,ie2 

where i and j represent single particle states of the first- and second-nuclei respectively and 
v is the effective nucleon-nucleon force. Formally (||) coincides with the corresponding term 
in the HF method. However in the DFM approach, the single-particle wave functions \i > 
and \j > are determined selfconsistently only inside each nucleus independently. To take into 
account the relative motion in the plane wave approximation one should multiply \i > by a 
factor exp(ikiVi) in the first nucleus and \ j > by exp(ik2Vj) in the second. Here k\ = K\/A\, 
&2 = K2/A2 and K a {a = 1,2) is the centre-of-mass moment of the first and second nucleus 
respectively. The relative momentum between a pair of nucleons of different nuclei is given by 
: K = k\ — &2- Its modulus can be defined globally as K 2 = 2mpE cm /h 2 or locally as ||: 
K 2 = 2mp(E — V — Vcoui)/fc 2 > where p = A\A2/(Ai + A2). In the local definition the set of 
equations: V = V(K); K = K(V) should be solved selfconsistently. 

The direct part of the HI potential is given by the first matrix element of @ and reads: 

V dir {D) = X d J df 1 df 2 p 1 (f 1 )p2(r2 - D)v(s), (4) 

where s = r\ — t%. The direct part of the HI potential eq.@ is just the convolution of the 
single-particle densities pi with the form factor v(s) of the central effective nucleon-nucleon 
force v. The constant X d is just the standard combination of the exchange parameters of the 
central nucleon-nucleon force: X d = w + b/2 — h/2 — m/4. This contribution is easily calculated 
numerically. It is important to note that the direct part of the HI potential || does not depend 
on energy if the nucleon-nucleon interaction v is energy independent. 

The finite-range exchange part of the HI potential is given by: 

V ex (D,K) =X e f dr 1 dr2p 1 (r 1 ,f 2 )p2(?2 - An - D)v{s)e iKs . (5) 

Here, as usual, X e = m + h/2 — 6/2 — w/4. To obtain this term a LDA for the density matrix is 
used. The simplest LDA is the SL approach, which in terms of centre-of-mass R = {r\ + 7*2) /2 
and relative s coordinates reads as: 

P(ri,r 2 ) = p{R) , (6) 
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where the Fermi momentum kp is related to the density by the usual Thomas-Fermi relation 
(|l]) (Here we use degeneracy factor 4 so the single particle density p is normalized to the mass 
number). This approximation is exact at HF level in infinite nuclear matter. To take into 
account finite size effects of the nuclei some modifications have to be made. Among different 
approaches to the DM, in the present paper we consider the Campi-Bouyssy (CB) @ one that 
is often used in DFM calculations. In the CB approximation the DM has the same form as the 
SL term (|3|), but instead of the Fermi momentum kp, uses an effective momentum k defined as: 

V = WT) iT{R) ~ i Ap(fl)) ' (7) 

In this equation r is just the exact quantal kinetic energy density. Consequently, to calculate it 
one needs to know the exact DM. To overcome this difficulty some approximate schemes for r 
are used, often they are based on its ETF value eq.(|2|) flTTf . 

Let us now derive the exchange term (||) using the ETF approach to the DM (8| for each 
isolated nucleus. This DM averaged over s can be written as: 

p{R,s) = ptf(R,s)+5 P (R,s), (8) 

where ptf is the Thomas- Fermi term, which is equivalent to the SL approximation (^). The 
averaged second-order semiclassical correction bp that we will use in this work reads as: 

Sp = —{Ap[j (k F s) [-xJoikps) - 3— }. (9 

72 kps p 6 kps 

In this approximation all the quantities are defined in a unique way. Eq.@ gives the ETF-?i 2 
correction to the DM in the case of a local HF potential. To take into account the non-locality 
of the single particle mean field, one should consider effective masscorrections in 5p. This case 
is analyzed in [|j in detail. It is also shown that the contribution of these non-local effects 
to the h 2 part of the DM are rather small as compared with the full h 2 correction if realistic 
finite-range forces are used. Therefore we will not consider the effective mass corrections to 
dp in this paper. We calculate the DFM potential using nucleon-nucleonforces with a Gauss, 
Yukawa and Coulomb-type form factors which cover most of the effective interactions used in 
DFM calculations. As typical examples we will consider the Gogny and M3Y forces ffl. 
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2.1 Gaussian-type force 

Let us consider a central nucleon-nucleon force with a Gaussian form factor: 

s 2 

v(s) = v exp( 2). (10) 

A* 

The exchange contribution to the DFM potential is given by: 

V ex (D) = X e v J d 3 r{Vsi(k[r}) + SV(k[r])}, (11) 
where the SL term Vsi is given by 

+ (<7 2 fcl + <7i/c 2 )(— + -3 ) ~ fc lMy + ^)]er/(-a) 
r , a 3 5a . . , , a 2 4 

+[ - fTia2( TV + ^ } + ( " 2 1 + ai 2)( V + V } 

-&iA;2^]-^Le;cp[— — a 2 ]}. (12) 



Here 0*1 2 is the sign of k\% (which are the Fermi momentum /cp of each nucleus related to the 
corresponding density through eq. ([!])), while a = K ± k\ ± &2. The summation in ( |i~2|) is taken 
over all possible a. The explicit form of (12) is given in [O]. To calculate the DFM potential in 



the CB approach we use the effective momentum (0) with the semiclassical ETF kinetic energy 
density (|). 

The second-order correction to the DFM is given by: 

OL 

Ufcl&2 ^(K + Xx) 2 /i 2 (JiT-Xi) 2 

+— T^K-ft + /2j(exp[ ] - exp[ ]) 



+(/i - /2)(exp[ 1 ] - exp[ ])] 

fel fc 2 2 V a 2 

1 /J 2 

H — 7=- Vo-icj 2 a"exp[ — -a 2 ]]}, (13) 



where a = .K" =F &i T &2> sci = &i + A; 2 , x 2 = fci — fc 2 and 

,(Vfc,) 2 _ ^ 2 4(V/^ 2 



h = A ki -2^ = ^[A Pi -;^] (14) 



ft = 2AA; l + — — = ^[aApi — ]. 

Let us now analyze these formulas. We will talk about the high energy limit when K/ki » 
l,i = 1,2 which is valid at high enough energy. However, it could also be valid at small energies 
in the outer region of nuclei when they overlap weakly. One can see that in this regime the SL 
term dominates and the h 2 correction is negligible. 

^■=- 5! Sf'' 3 ^ fee "^ +o «|i 2 i)- < i5 > 

On the contrary, at small energies the h 2 correction plays an important role. The SL term in 
the low energy limit, e.g. when ¥- << 1, looks like 



+— [*(exp[-£(ki - k 2 ) 2 ] - exp[-^(h + k 2 ) 2 ]) 



+Uh + hfexpK-^ih + k 2 ) 2 } - - k 2 fexp[{-^-{k 1 - k 2 f] 



\k x k 2 {exp\-^{k x - k 2 f\ + expl-^ih + k 2 ) 2 })}} + 0([^] 2 ). (16) 



This formula was obtained in [13| for the static SL DFM potential with the Gogny force. The 
low energy limit of the h 2 correction is: 

SV(R) = {[3(| + fml + k 2 )-±]- 2(k lfl + k 2 f 2 )) 

x [ exp [-VL( kl + k2 f] _ expl-^ih - k 2 ) 2 }} 
-3(g + + k 2 f exp[-^(fc! + k 2 ) 2 ] 

-(k.-k^eM-^ih-^) 2 }] 
-^hfaKh + k 2 ){h + f 2 )exp[—^-(ki + k 2 f] 

+(h - k 2 )(h - f^expi-^ih - k 2 ) 2 ]]} + 0((^) 2 )- (17) 

In a realistic calculation both limiting cases: K/ki << 1 and ki/K « 1 could take place at 
the same relative energy, but at different points in space. 
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2.2 Yukawa-type force 



For a Yukawa-type force 



e -(3s 



«(*) = "o^- (18) 



one will get 



_ a 2 /? 3 a 4 /3 /? 5 . 

+ L[-^(^--^ r - M ) 

. , , w a/3 3 a 3 /3. , , ,a 2 3 /3 s ,. .a. 
+(cr 2 A;i + aik 2 ){— — ) + hk 2 {— — )\arctg{-) 

. . a 5 a(3 A a 3 2 . 

+ [ - aia2{ 120 + -24 ~ 

a 4 a 2 /? 2 B 4 aB 2 a 3 

+(a 2 h + a^X^ - - JL + y + hk 2 (J±- - -M(a)}, (19) 



where only the combination 

f ds 1 a 2 + j3 2 

Jf (ai) — J\(oi 2 ) = / — \cos{a\s) — cos(a 2 s)]exp(—3s) = -In—* 

J s 2 af + p z 

appears. Again if one uses the CB effective momenta one will get the CB approximation to the 
DFM. The explicit expression of ( |l9| ) can be found in []l2j ]. The second-order correction to the 
DFM potential is: 

+ hh) - + f (k 2 + k 2 ) - ^ 2 )(fi + g)] £ *i<*«rc^) 

+ [f (*l/l + fc*0 + - ^(*? + ^2 + /3 2 ))(| + g) 

(^ + x 1 ) 2 + /3 2 , (K- Xl ) 2 + B 2 ^ 
X[Ln (K + x 2 ) 2 + 3 2+ (K-x 2 ) 2 + 3 21 

+ ^A + ^)-^ + *S)(g + g)]- §±|^g 

Considering the high energy limit, where the Ti 2 correction to the potential is negligible, one will 
get 

4-7T -> h ■ 

Vsi = -^Pi(r)p 2 (r - R) + 0([-^] 4 ). (22) 
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These results for a Yukawa force correspond to the convolution of the densities of the isolated 
nuclei with a zero-range pseudopotential and an extra 1/E energy-dependence. Therefore at high 
enough energy, when K/ki > 1, a simple zero-range pseudopotential can be used for calculating 
the exchange part of the DFM potential in this case. The parameters of this pseudopotential 
are determined by the range (3 of the effective nucleon-nucleon force (here M3Y). The inverse 



energy dependence in (22) should be taken into account. 

Now we consider the low energy limit where h 2 correction is more important. The SL term 
in this case looks like: 

+ ^i^kfkl + - k 2 ) 2 - g + - klfikl + k\ - p 2 )]] 
- f P ? + kl)arctg(^±^) - (kf - kl)arctgf^)] 
+ (^ + ^(k 2 + k 2 ) V ^?)ln {h + k2)2 + (32 ] I C(( K ?) 
The corresponding fi 2 correction reads: 

™ = -9^^ fe <| + g' (24) 
2A-i fro 

- ' -[/3 2 (fri/i + fca/a) + (k\ - kDikxh ~ k 2 h)\ 



Q 

1 \fu f i i n 3 ^2 , t.2 , o2\,9l . 92^ 17 ( fc i + k 2 ) 2 + f3 2 K. 
2 Pi/i + W " j (*! + ^2 + /? )(^ + T 2 )]ln ( kl - k 2 ) 2 + P } + °%- 



where 



Q = /3 4 + 2/^(fc 2 + fr^) + (fc 2 -fr^ (25) 



2.3 Coulomb force 

Let us now consider the Coulomb force: 



v(s) = £ (26) 



In this particular case one will get the explicit expression for the SL term: 

^ e2 ^ 2 {lif 3 fc 1 fc2 + -( 

Kk\kl 1 30 30 



Vsi = - 9 -^^{^k lk2 + ^( k 2 + k 2 )Kk lk2 
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+ [ 240 24 {kl + k2 > 16 (fcl k2) l[Ln (K + h- k 2 f +ln (K-k 1 + k 2 )^ 
-\^{k\ + kl) + ^(h + k 2 f[( kl - k 2 ) 2 - hk 2 ]]ln^±^±^ 

+ [^(kf ~ k 3 2 ) + ±(h - k 2 f[(h + k 2 ) 2 + kl k 2 ]]ln ^±^^ }, (27) 

where now pi and k{ (i=l,2) correspond to the density and Fermi momentum of protons in 
each nucleus: pi = k{ = ki^ p . On first sight of ( p5| ) one can think that at K = ±/ci ± k 2 
this formula contains divergent terms. However, one can see that the corresponding terms that 
appears twice cancel each other. 

The second-order correction to the DFM potential is: 

+[§ (fci/i + k 2h) + {^r-f(k 2 + fc 2 2 ))(| + 1) 

+ lkkfh + k 2 f 2 ) -h$ + + ^)]/n ( ^ + Xl)2 



2 v-i,i ■ -w, 2 v»r ■ 1 k 2 n "\K-x x f 

5«A-«A)-i«-«)(| + |)li-^. 



Considering the case X = a?i 2 one can see that ( ^8|) still remains well defined. In the high 
energy limit we get again the simple formula: 

%i = -^Pl(r> 2 (r-l?) + 0([§] 4 ) (29) 
In the low-energy limit one obtains 

= -^^^ + ^ " ^ " k2)2ln ^W } + 0(( f )2) (30) 

and its fi 2 correction reads as: 

In the case fci = &2 one should take care about the correct consideration of the limit k\ — ► k 2 ■ 
In this case one should start from eq.(26) for the Yukawa force, put k\ = k 2 and make the limit 
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Plimn-^o- After path integration one will get at k\ = k 2 : 

SV = -^[(Vhf + (Vk 2 n (32) 

It is very important to point out that the SL, CB or ETF approximations to the DM described 
in this section are used for calculating the exchange part of the DFM potential irrespective of the 
way that the density of each isolated nucleus is obtained. In the following we will refer to these 
potentials as SL, CB or ETF approaches to the DFM potential. Notice that for determining 
completely the DFM potential, the density of each isolated nucleus that enters through k\ and k 2 
in the previous formulas as well as the effective nucleon-nucleon force used have to be specified. 



3 Results 

In this Section we present the results for the DFM potential obtained using the different prescrip- 
tions considered in Section 2. The nuclear part of these potentials calculated with the realistic 
effective nucleon-nucleon force analyzed in this paper is attractive. Hence the relative momen- 
tum K in the local definition is always real when the energy exceeds the Coulomb barrier. On 
the other hand, the "frozen density" approximation, used in DFM, is still valid for fast enough 
collisions or, in other words, at rather high energies. The Coulomb barrier for the ie O — 16 O 
system is about lOMeV. Therefore, we will consider only collisions above the Coulomb barrier 
where K is always real. 

To check the validity of the ETF DM to the calculation of an ion-ion potential, we consider 
the exactly soluble HO model for the isolated nuclei. This model consists of describing the 
ground state of each colliding nucleus using HO wavefunctions. In this case the exchange part 



of the potential is easily calculated due to factorization of variables in the DM [11]. 

First of all, let us briefly describe the effective nucleon-nucleon force used in this paper. Here 
we will consider the Gogny force that is used in nuclear structure calculations and a MY3-type 
force that is typical in DFM calculations. 

The effective Gogny force [|l0| consists of a finite-range Brink-Boeker term together with a 
density-dependent zero-range term 

..2 

t4 



2 2 

v (r) = J2[ Wi + hP° - h { P T - mi p° p-] eX p{- r —) + t 3 {l + x 3 P <T )p 1/3 ( !1 ^)5(r), (33) 
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where P a and P T are the usual spin and isospin exchange operators, the coefficients Wi,bi,hi 
and mi are the parameters of the central force and pi the range of the Gaussian form factor. The 
values of these parameters are given in |T(J. In DFM calculations the density dependent term is 
taken in the so-called "sudden" approximation, which implies that the density of the composite 



system is simply the sum of individual densities: p(r) = pi(r) + p2(r — R) [13|. 
The density-dependent M3Y force is defined as 

v(r)=F(p)g(E)v M 3Y(r), (34) 

where VM3Y( r ) is the finite range of the M3Y Paris j^] or Reid-Elliott [p^f] force with a Yukawa- 
type formfactor. For a BDM3Y force the density-dependent term is given by j|: 

F BD {p) = C{l-ap 1 ), (35) 

while the DDM3Y force 0] reads: 

F DD = C(l - ae-T"). (36) 

Again, the density-dependent term is considered in the "sudden" approximation. The energy 
dependent term is taken as g(E) = 1 — CeE. 

In this paper we are going to consider the validity of different approximations to the DM for 
the calculation of the DFM potential rather than obtaining a realistic potential for describing 
the elastic scattering data. In order to do this wee will make the two following simplifications. 

The first concerns the way that the ground state density of each isolated ie O nucleus is cal- 
culated. To do this, we minimize the quantal ground state energy with respect to the oscillator 



parameter e = y/mojfh. In the semiclassical calculations using the ETF, SL and CB approaches 



to the DM we first write the ground state energy using the Kohn-Sham [16| scheme (see refer- 
ence H for more details) and then minimize with respect to the HO parameter. The binding 
energy and root mean square radius (RMSR) of the ls O ground state obtained quantally and 
semiclassically (with the ETF, SL and CB approaches) are collected in Table 1. One can see 
that in all cases the results obtained with the ETF approach to the DM give a better descrip- 
tion of the quantal ground state properties of 16 O than that obtained using the semiclassical SL 
or CB approximations to the DM. On the other hand, among the considered M3Y forces, the 
BDM3Y1* (Paris) force gives a RMSR that is very close to that obtained with the Gogny force 
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and consequently the same ground state density for the ie O. Therefore, we will use the Gogny 
and BDM3Y1* (Paris) forces in our analysis. 

The second remark deals with the definition of relative momentum K. If one uses the local 
definition of K and solves the system of equations: V=V(K); K=K(V) selfconsistently, the 
nuclear and Coulomb potentials are not separated. In this paper we are going to consider the 
different kinds of forces independently. For this reason and for the sake of simplicity we will use 
the global definition of K. 

Now we consider the resulting ie O — 16 O ETF potentials calculated using the quantal HO 
density for each ie O nucleus and the Gogny and BDM3Y1* forces. In Figure 1 these potentials 
calculated at several energies are displayed by dashed and solid lines respectively. The upper 
dashed line is the energy independent sum of the Brink-Boeker and zero-range direct parts of 
the DFM potential in the Gogny force case. Thus, the difference between the full potential and 
this energy independent direct part gives the contribution of the one-particle exchange due to 
the finite (energy dependent) and zero (energy independent) ranges of the Gogny force. One can 
see that at low energies the total potential is purely attractive due to the exchange interaction 
coming from the finite-range part of the Gogny force. When the energy increases, the total 
DF potential tends to its energy independent part (direct Brink-Boeker and full zero-range 
contributions). This latter potential is repulsive at small distances and attractive when the 
separation between the nuclei is increased and tends to zero at large distances. 

As can be seen from eq.([34|), the BDM3Y1 force depends linearly on the energy. This force 
has two essential differences when compared with the Gogny interaction. First of all, in this 
case the direct DFM potential (|3|) is energy dependent. The upper solid line corresponds to this 
direct potential in the E = case. On the other hand, taking into account the linear energy 
dependence of the BDM3Y force and the inverse energy dependence of the exchange potential 



(22) (in the high-energy limit) one sees that the exchange potential in this case does not vanish 
when the energy tends to infinity and, consequently, the total potential does not tend to the 
direct part in this limit. The total and direct ETF DFM potentials obtained with the BDM3Y1* 
force are also displayed in Figure 1. Looking at this figure one can see that in the region of 
energies of physical interest (E cm < lOOOMey) the total DFM potentials obtained using quantal 
HO densities and Gogny and BDM3Y1* forces are very similar. At higher energies they start 
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to show differences. However, these differences should be considered only in the mathematical 
aspect because the nucleus-nucleus potentials at these high energies can not be defined correctly. 

Now we turn to the discussion of the validity of different approximations to the DM for 
calculating DFM potentials considered in the present paper. As was pointed out above, the ex- 
change effects are larger at small energies and their contribution falls when the energy increases. 
Therefore it is reasonable to expect that the influence of the different approaches to DM on the 
DFM potential will be more important at low energies. Hence to investigate the applicability 
of the different approximations to the DM for calculating the DFM potential, we will consider 
the static case (K=0). To this aim we define a function cr(D), which is the deviation of the 
approximate total DFM potential with respect to the exact potential: 

a ^ = VSL,CB,ETF{D) - Vexact(D) ^ 

where V exac t{D) is the DFM potential calculated with the quantal HO DM at a separation 
distance D, whereas Vsl, Vcb or Vetf are the DFM potentials calculated with the SL, CB or 
ETF approximations. We perform two sets of calculation using two different densities for the 
isolated ie O nucleus to obtain the Vsl,Vcb and Vetf for the Gogny and BDM3Yl*forces. 

In the first case we used the densities that minimize the ground state Kohn-Sham energy 
calculated in the SL, CB or ETF approach for an uncharged ie O nucleus. In other words this 
case corresponds to the selfconsistent calculation, where the same approach to the DM is used 
to obtain the DFM potential and the density of each isolated nucleus. 

These results are plotted in Figures 2 and 3 for the Gogny and BDM3Y1* forces respectively. 
These figures clearly show that the DFM potential Vetf calculated using the selfconsistent ETF 
densities reproduces better the quantal DFM potential than Vsl or Vcb calculated with the 
selfconsistent SL or CB densities. This result is obvious if one takes into account that for an 
isolated 16 O nucleus the selfconsistent density in the ETF approximation is closer to the quantal 
density as compared with the other approximation to the DM, as can be seen from Table 1. 

In the second set we obtain the different DFM potentials using the quantal HO density for 
each ie O nucleus and the Gogny and BDM3Y1* forces. The results for the function cr(D) are 
plotted in Figures 4 , 5 and 6 for the Gogny, BDM3Y1* and Coulomb forces respectively. 

In this case the situation is more complicated. From Figures 4 and 5 one can see that Vsl 
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(dotted line) does not reproduce the quantal potential for distances D > 3fm. The potential 
Vetf (dashed-dotted line) is very close to the quantal DFM potential for D < 6fm. If D > 6fm, 
Vcb (dashed line) gives the best agreement with the quantal DFM potential, while for D < 6fm 
it works slightly worse than Vetf- At large distances R «10 fm, the nuclear part of the nucleus- 
nucleus potential is very small wh compared with its Coulomb part. Consequently, the rather 
strong deviation ( within 10%) shown by the the nuclear part is not very important. However, 
from the formal point of view one should take it into account. 

As is discussed in M all the considered approximations to the DM (SL, ETF or CB) are 
actually distributions. They are very efficient for obtaining expectation values, but there is no 
reason for reproducing the quantal DM at each point in coordinate space. If one looks for the 
asymptotic behaviour of the different approximations to the DM considered in this paper, one 
can see that at finite s, the SL and CB DM behave as : p(R,s) oc 2Ca 2 R 2 exp[— a 2 R 2 ] when 
R — > oo (more precisely when aR >> 1,R >> s), where C is the normalization constant. The 
ETF DM has a different asymptotic behaviour: p(R,s) oc 2Ca 2 R 2 (l - ^a 4 s 2 R 2 )exp[-a 2 R 2 ]. 
At the same time the quantal DM asymptotically behaves as p(R,s) oc 2Ca 2 R 2 exp[— a 2 R 2 ]. 
One can see that at very small s and large R the SL, CB and ETF DM are similar to the 
quantal density matrix. However, when s is increased, the ETF DM starts to deviate from the 
quantal values. Due to the extra R 2 dependence, the ETF DM gives a worse approximation to 
the quantal DM as compared with SL and CB DM if one uses the same quantal HO density for 
obtaining the different approaches to the quantal DM. 

This lack of accuracy of the ETF DM in describin the quantal DM is not important when 
one calculates the ground state properties of isolated nuclei. In the region of wrong asymptotic 
behaviou the DM is too small to give any significant contribution to the total energy as can be 
seen from Table 1. However, for calculating DFM potentials the situation is different. The values 
of densities over all the nucleus contribute to the DFM potential at large distances Therefore, the 
wrong asymptotic part of the ETF DM can give a significant contribution to the DFM potential 
at large separaon distances. Thus, we come to the conclusion that the DFM potential is a more 
sensitive tool for investigating the one-particle density matrix in nuclei than fully integrated 
quantities such as the binding energy or RMSR. 

For improving Vetf calculated with the quantal HO density at large separation distances 
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D, we can make the following modification of the ETF DM. We introduce a cutoof radius R co 
that removes the h 2 correction from the SL DM. Therefore this new ETFM DM is is just the 
ETF DM when R < R co and the SL DM if R > R co and has the right asymptotic behaviour 
for large R values. The results of the DFM potential calculated with this new ETFM DM are 
displayed in Figures 4, 5 and 6 by solid lines. It is clear that this new approximation gives a 
better description of the exact DFM potential at all distances D. However, a new parameter, 
the cutoof radius in the outer region, is introduced in this case. The solid line corresponds to 
the approximated DFM potential obtained with this modified ETF DM (Vetfm), where the h 2 
contribution has been cut off after R co . One can see that Vetfm agrees to within 1% with the 
quantal DFM potential. 

In Figure 6 we show the results for the long-range Coulomb force. For D < 4fm Vsl gives a 
better result than Vcb and at D > 4fm it starts to deviate from the t quantal DFM potential 
values. The Vetf gives a very good description of the quantal DFM potential and gives an 
error within 0.25% if 6 < D < 12fm. Furthermore, for the Coulomb case the use of Vetfm 
significantly improves the agreement with the quantal DFM potential in all of the considered 
region. 

It is also interesting to consider in our model calculations the height and position of the 
Coulomb barrier. The results corresponding to the selfconsistent calculations are collected in 
Table 2. In this case the ETF DFM potential gives a clearly better approximation to the exact 
quantal values when compared with SL and CB DFM approaches. This result is obvious if one 
takes into accountthe fact that the selfconsistent density for the ETF case is very close to the 
quantal density. The height and position of the Coulomb barrier calculated with the different 
approximations to the DFM potential obtained with the same quantal density are shown in 
Table 3. In this case one can see that the better agreement with the quantal results is obtained 
using the CB approach to the DFM potential. To improve the ETF DM potential we again use 
the ETFM approxiationto the DM where a cutoof radius R co is introduced. First of all,from 
Fig.4-6 one can obtain the following conclusion: the longer the range of the force, the smaller the 
value of cutoof radius that should be used to reproduce the quantal DFM. The Coulomb barrier 
comes from the superposition of longe-range Coulomb force (R co = 5/m) and a shorter-range 
nuclear force. We choose the cutoof radius R co = 5. 5/m and R co = 5/m for the Gogny and 
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BDM3Y1* forces respectively. One can see that the ETFM DFM potential nicely reproduces the 
quantal result for both forces. From Table 3 one can also see thatthe the Gogny and BDM3Y1* 
interactions give very close results for the Coulomb barriers. 

4 Summary 

In this paper we have derived the one-particle exchange contribution to the DFM potential 
using the recently proposed ETF approach to the DM. We consider here nucleon-nucleon forces 
with the widely used Gauss, Yukawa and Coulomb form factors. As representative examples of 
these forces, we have performed our numerical analysis with the Gogny and BDM3Y1* effective 
interactions. 

To obtain the densities of each isolated 16 nucleus, we have used trial HO wavefunctions. 
The oscillator parameter e = ^muj/h is obtained by minimizing the exact ground state energy 
in the quantal case and the corresponding Kohn-Sham energy when the quantal DM is replaced 
by its SL, CB or ETF approximations. 

If each self-consistent density (SL, CB or ETF) is used for obtaining its corresponding DFM 
potential, it is found that the ETF potential gives a better description of the quantal one than 
the SL and CB DFM potentials. 

However, if one uses the same selfconsistent quantal density to obtain all the approximated 
DFM potentials, the situation is more complicated. First of all, it is clearly seen that the SL 
approach to the DFM potential compares worse with the quantal potential than the CB or 
ETF approaches. In all the analyzed cases the ETF approach to the DFM potential gives 
a better description of the quantal potential at small and intermediate separation distances 
(D < 6fm) when compared with the potential obtained in the CB approximation. However, at 
large separation distances, the ETF approach fails in describing the quantal potential. From a 
physical point of view, this is not very important because the nuclear part of the total ion-ion 
potential is very small at these distances. The reason for this deviation of the ETF potential 
from the quantal potential is due to the different asymptotic behaviour of the quantal and ETF 
DM of each isolated nucleus. To improve the ETF DM of each nucleus at large distances, we 
drop its ^-contribution after a cutoof radius R co - In this way one obtains a modified DM which 
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coincides with the full ETF DM at R < R co and with its SL part DM if R > R co . We have 
numerically shown that this ETFM DM produces a DFM potential in very good agreement with 
the quantal result in the whole range of distances. 

Hence we reach the conclusion that the DFM potential is a sensitive tool for investigating the 
quality of the approximations to the DM in the whole range of R- values. While in the calculations 
of ground state properties of nuclei the asymptotic behaviour of the DM is suppressed by the 
small values of the particle density at large R, in the calculation of the DFM potentials it is 
clearly revealed. 
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6 Table Captions 



Table 1. Binding energy and root mean square radius < r 2 y- 1 / 2 of an uncharged 16 nucleus 
calculated with trial HO wavefunctions quantally (QM) and using the SL, CB and ETF approxi- 
mations to the density matrix. The results obtained with the Gogny and several M3Y forces are 
displayed. The centre-of-mass correction to the energy has been taken into account. The HO 
parameter which minimizes the binding energy for each force is given by a = 3/2 < r 2 >~ 1 / 2 . 

Table 2. Height (Vb) and position (Rb) of the Coulomb barrier calculated with the quantal 
DFM potential (QM) and with its SL, CB and ETF approximations using the corresponding 
selfconsistent densities for the Gogny and BDM3Y1* forces. 

Table 3. The same as Table 2 but with the SL, CB and ETF DFM potentials obtained with 
the quantal density. 



19 



Figure 1. Direct and exchange contributions at several energies to the ETF DFM potential 
calculated with the quantal density. Solid and dashed lines correspond to the results obtained 
with the BDM3Y1* and Gogny forces respectively. 

Figure 2. Relative deviation from the quantal DFM potential with respect to its SL (dotted), 
CB (dashed) and ETF (solid) approximations obtained with the corresponding selfconsistent 
density and the Gogny force. 

Figure 3. The same as Figure 2, but for the BDM3Y1* force. 

Figure 4. Relative deviation from the quantal DFM potential with respect to its SL (dotted), 
CB (dashed) , ETF (dashed-dotted) and ETFM (solid) approximations using the quantal density 
and the Gogny force. 

Figure 5. The same as Figure 4, but for the BDM3Y1* force. 

Figure 6. The same as Figure 4, but for the Coulomb force. 
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Table 1 
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< r 2 >]il 
^ ' ^CB 


P i T7"~T~< T7< 


< r 2 >!/' 
^ ' ^ ETF 


Gogny 
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2.63 


-129.49 


2.67 


-133.80 


2.65 


-134.35 


2.66 


BDM3Y1 
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2.55 


-112.67 


2.57 


-114.63 


2.57 


-115.76 


2.56 


BDM3Y2 
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2.55 
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BDM3Y3 
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BDM3Y1* 
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2.63 


DDM3Y* 
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-117.13 
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Table 2 





QM 


SL 


CB 


ETF 


Gogny 


R B (fm) 


8.36 


8.51 


8.42 


8.38 


Y B {MeV) 


10.35 


10.18 


10.27 


10.35 


BDM3Y1* 


Rb 


8.46 


8.57 


8.50 


8.43 


V b 


10.18 


10.06 


10.12 


10.24 



Table 2 





QM 


SL 


CB 


ETF 


ETFM(R co =5fm) 


Gogny 


R B (fm) 


8.36 


8.40 


8.35 


8.35 


8.37 


Y B (MeV) 


10.35 


10.31 


10.36 


10.38 


10.35 


BDM3Y1* 


R B (fm) 


8.46 


8.51 


8.45 


8.40 


8.45 


Y B (MeV) 


10.18 


10.13 


10.19 


10.25 


10.19 
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